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Abstract — This paper presents the results of a nonlinear Finite Element (FE) analysis conducted on Reinforced High 
Strength Concrete (HSC) beams strengthened with Glass fibre reinforced polymer (GFRP) laminates. Modeling the 
complex behavior of reinforced concrete, which is both non-homogeneous and anisotropic, is a difficult confront in the 
finite element analysis of civil engineering structures. The accuracy and convergence of the solution depends on factors 
such as mesh density, constitutive properties of concrete, convergence criteria and tolerance values etc., Three- 
dimensional finite element models were developed using a smeared crack approach for concrete and three dimensional 
layered elements for the FRP composites. The results obtained through finite element analysis show reasonable 
agreement with the test results. 

Keywords — Reinforced Concrete, Nonlinear Analysis, Finite Element Analysis and GFRP. 

I. INTRODUCTION 

Strengthening or upgrading becomes necessary when the structural elements cease to provide satisfactory strength 
and serviceability. Fiber Reinforced Polymer (FRP) composites can be effectively used as an external reinforcement for 
upgrading such structurally deficient reinforced concrete structures. The most common types of FRP are aramid, glass, and 
carbon; AFRP, GFRP, and CFRP respectively. Many researchers have found that FRP composites applied to reinforced 
concrete members provide efficiency, reliability and cost effectiveness in upgradation[4,5,9]. The use of FEA has been the 
preferred method to study the behaviour of concrete. Wolanski[2] studied the flexural beha\ ior of reinforced and prestressed 
concrete beams using finite element analysis. The simulation work contains areas of study such as Behavior at First 
Cracking. Behavioral Initial Cracking. Behavior beyond First Cracking. Behavior of Reinforcement Yielding and Beyond, 
Strength Limit State, Load-Deformation Response of control beam and Application of Effective t i It in t 

Deflection, Decompression, Initial Clacking. Secondary Linear Region. Behavior of Steel Yielding and Beyond, Flexural 
Limit State of prestressed concrete beam. Arduini, et al.|3] used finite element method to simulate the behaviour and failure 
mechanisms of RC beams strengthened with FRP plates. The FRP plates were modeled using two dimensional plate 
elements. However the crack patterns were not predicted in that study. Kachlakev, etal.[7] studied the finite element 
modeling of reinforced concrete structures strengthened with FRP Laminates with ANSYS and the objective of this 
simulation was to examine the structural behaviour of Horsetail creek bridge with and without FRP laminates and establish a 
methodology for applying computer modeling to reinforced concrete beams and bridges strengthened with FRP laminates. 

II. EXPERIMENTAL BEAMS 

2.1Materials used 

The concrete used for all beam specimens had a compressive strength of 64MPa. The concrete consisted of 450 
kg/m3 of ordinary Portland cement, 780 kg/m3 of fine aggregate, 680 kg/m3 of coarse aggregate, 450 kg/m3 of medium 
aggregate, 0.36 water/cement ratio and 0.8% of hyperplasticizer. The longitudinal reinforcement consisted of high yield 
strength deformed bars of characti u tii tri n th ' • i II he lal ral ties consisted of mild steel bars of yield strength 
300MPa. The specimens were provided with 8mm diameter stirrups at 150 mm spacing. Two types of GFRP laminates were 
used for the study, namely, Chopped Strand Mat (CSM) and Uni-directional Cloth (UDC) of 3mm and 5mm thickness. The 
properties of GFRP are shown in Table 1. 

TABLE1 Properties of GFRP Laminates 



Type of 
GFRP 


Thickness 
(mm) 


Elasticity Modulus 
(MPa) 


Ultimate Elongation (%) 


Tensile Strength 
(MPa) 


Chopped Strand 
Mat 


3 


7467.46 


1.69 


126.20 


5 


11386.86 


1.37 


156.00 


Uni-Directional 
Cloth 


3 


13965.63 


3.02 


446.90 


5 


17365.38 


2.60 


451.50 
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2.2 Details of Beams 

A total of 15 beams were tested. The main test variables considered in the study were steel reinforcement ratio, 
type of GFRP laminate and thickness of GFRP laminate. The beams were 150 x 250 mm in cross-section and 3000 mm in 
length as shown in Figs. 1-3. The beams of A series were reinforced with two numbers of 10 mm diameter bars giving a steel 
ratio of 0.419%. The beams of B series were reinforced with three 10 mm diameter bars giving a steel ratio of 0.628%. The 
beams of C series were reinforced with three 12 mm diameter bars giving a steel ratio of 0.905%. Stirrups of 8 mm diameter, 
at a spacing of 150 mm, were used for the beams. Out of fifteen beams, three served as control beams and the remaining 
beams were strengthened with GFRP laminate. The details of beams are presented in Table. 2 




Fig.2 Details of 'B ' Series Beams Fi S 3 Details of *C Series Beams 



TABLE 2 Specimen Details 

GFRP Laminate 



RA 


0.419 




RAC3 


0.419 


CSM 


RAC5 


0.419 


CSM 


RAU3 


0.419 


UDC 


RAU5 


0.419 


UDC 


RB 


0.628 




RBC3 


0.628 


CSM 


RBC5 


0.628 


CSM 


RBU3 


0.628 


UDC 


RBU5 


0.628 


UDC 


RC 


0.905 




RCC3 


0.905 


CSM 


RCC5 


0.905 


CSM 


RCU3 


0.905 


UDC 


RCU5 


0.905 


UDC 



3.183 
1.909 
3.183 

1.909 
3.183 
1.909 
3.183 



2.3 GFRP Laminate Bonding Technique 

Glass Fibre Reinforced Polymer (GFRP) laminates were used for strengthening the beams. The soffit of the beam 
was well cleaned with a wire brush and roughened with a surface grinding machine. Two part epoxy adhesive consisting of 
epoxy resin and silica filler was used to bond the GFRP laminates. The adhesive was spread over the beam soffit with the 
help of a spread. The GFRP laminate was applied gently by pressing the sheet from one end of the beam to the other along 
the length of beam. 

2.4 Experimental Test Set-up 

All the beams were tested under four point bending in a loading frame of 750 kN capacity. The effective span of 
the beam was 2800 mm with 100 mm bearing at the ends. The deflections were measured at mid-span and load-points using 
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dial gauges of 0.01 mm accuracy. The crack widths of beams were measured using a crack detection microscope with a least 
count of 0.02 mm. Figure. 4 shows the loading arrangement and instrumentation adopted for the test. 



13=VHiH 





Fig.4 Experimental Test Set-up 

III. FINITE ELEMENT MODELING 

3.1 Element types 

3.1.1 Reinforced Concrete 

Solid65 element was used to model the concrete. This elemenl has eight nodes with three degrees of freedom at 
each node - translations in the nodal x, y, and z directions. This element is capable of plastic deformation, cracking in three 
orthogonal directions, and crushing. A schematic of the element is shown in Figure 5 [6] . Smeared cracking approach has 
been used in modeling the concrete in the present study [81 

3.1.2 Reinforcement 

The geometry and node locations for Link 8 element used to model the steel reinforcement are shown in Figure 6. 
Two nodes are required for this element. Each node has three degrees of freedom, translations in the nodal x, y, and z 
directions. The element is also capable of plastic deformation. 

3.1.3 FRP Composites 

A layered solid element, Solid46, was used to model the FRP composites. The element has three degrees of 
freedom at each node and translations in the nodal x, y, and z directions. The geometry, node locations, and the coordinate 
sy stem are shown in Figure 7. 
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Fig.6 Link 8 - 3-D spar element 
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Fig. 7 Solid46 - 3-D layered structural solid 



a for material properties as follows: 



e,/r)Mpa 



te, ANSYS requires input dai 

• Elastic modulus (Ec ) MPa 

• Ultimate uni-axial compressive strength (f'c) Mpa 

• Ultimate uni-axial tensile strength (modulus of ruptur 

• Poisson's ratio (D =0.2) 

• Shear transfer coefficient (fit) which is represents conditions of the crack face. The value of/?? ranges from 0.0 tc 
1.0, with 0.0 representing a smooth crack (complete loss of shear transfer) and 1.0 representing a rough crack (nc 
loss of shear transfer) [I] . The shear transfer coefficient used in present study varied between 0.3 and 0.4 

• Compressive uni-axial stress-strain relationship for concrete 

For steel reinforcement stress-strain curve for the finite element model was based on the actual stress-stain curve obtained 

from tensile test. 

■Material properties for the steel reinforcement are as follows: 

• Elastic modulus (Es) Mpa 

• Yield stress(/y) Mpa 

• Poisson's ratio (□) 

Material properties for the GFRP laminates are as follows: 

• Elastic modulus 

• Shear modulus 

• Major Poisson's ratio 
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Fig.8 Simplified Uniaxial Stress-Strain Curve for Concrete Fig.9 Stress-Strain Curve for Steel 

Reinforcement 
3.2 Cracking of Concrete 

The tension failure of concrete is characterized by a gradual growth of cracks, which join together and finally 
disconnect larger parts of the structure. It is a usual assumption that crack formation is a brittle process and the strength in 
the tension-loading direction abruptly goes to zero after big cracks or it can be simulated with gradually decreasing strength. 
The cracked concrete material is generally modeled by a linear elastic fracture relationship. Two fracture criteria are 
commonly used, the maximum principal stress and the maximum principal strain criterions. When a principal stress or strain 
exceeds its limiting value, a crack is assumed to occur in a plane normal to the direction of the principal stress or strain. Then 
this crack direction is fixed in the next loading sequences. In this study the smeared-crack model was used. A three- 
dimensional failure surface for concrete is shown in Figure 9. 
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Fig.9 3-D Failure surface for 
3.3 Finite element discretization 

The finite element analysis requires meshing of the model. For which, the model is divided into a number of small 
elements, and after loading, stress and strain are calculated at integration points of these small elements. An important step in 
finite element modeling is the selection of the mesh density. A convergence of results is obtained when an adequate number 
of elements is used in a model. This is practically achieved when an increase in the mesh density has a negligible effect on 
the results. 



3.4 Non-linear solution 

In nonlinear analysis, the total load applied to a finite element model is divided into a series of load ii 
called load steps. At the completion of each incremental solution, the stiffness matrix of the model is adjusted to reflect 
nonlinear changes in structural stiffness before proceeding to the next load increment. The ANSYS program (ANSYS 2010) 
uses Newton-Raphson equilibrium iterations for updating the model stiffness. Newton-Raphson equilibrium iterations 
provide convergence at the end of each load increment within tolerance limits. In this study, for the reinforced concrete solid 
elements, convergence criteria were based on force and displacement, and the convergence tolerance limits were initially 
selected by the ANSYS program. It was found that convergence of solutions for the models was difficult to achieve due to 
the nonlinear behavior of reinforced concrete. Therefore, the convergence tolerance limits were increased tc 
times the default tolerance limits in order to obtain convergence of the solutions. 

TABLE 3. Material Models for the Calibration 



Element Type 



Material Properties 



Linear Isotropic 



Multi Linear Isotropic 




Stress 


Strain 


Point 1 


19.2 


0.00048 


Point 2 


39.68 


0.0019 


Point 3 


54.01 


0.0027 


Point 4 


58.56 


0.0029 


Point 5 


64 


0.0032 



ShfCf-Op 



UnCompSt 
BiCompSt 



BiCompSt 

UnTensSt 



Linear Isotropic 

I 7467.46Mpa ~ 
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Fig.16. Flexural Crack Signs 



Fig. 17. Concrete Crack Signs 



IV. RESULTS AND DISCUSSION 

it results at first crack, yield and ultimate stage of non-strengthened and strengthened beams. 
TABLE 4 Principal Results of Tested Beams 



Beam Designation 


Loading Stages of Beams 


First Crack Stage 


Yield Stage 


Ultimate Stage 


Per (kN) 


□ CT (mm) 


Py(kN) 


D y (mm) 


P„ (kN) 


, u (mm) 


RA 


14.39 


1.26 


29.42 


7.91 


41.68 


21.05 


RAC3 


16.52 


1.41 


36.77 


9.02 


51.48 


33.46 


RAC5 


21.28 


3.67 


46.58 


10.1 


66.19 


46.81 


RAU3 


32.94 


7.98 


51.48 


11.42 


71.09 


53.26 


RAU5 


36.81 


9.23 


53.7 


10.74 


78.45 


57.21 


RB 


28.32 


3.68 


39.22 


8.11 


53.93 


31.28 


RBC3 


30.95 


4.71 


51.48 


11.35 


61.29 


36.23 


RBC5 


32.17 


4.97 


53.24 


12.41 


63.74 


56.91 


RBU3 


33.69 


9.35 


58.8 


12.85 


88.25 


61.04 


RBU5 


39.41 


11.14 


63 


12.69 


100.51 


65.59 


RC 


30.37 


4.45 


44.13 


10.19 


58.84 


37.89 


RCC3 


31.88 


3.85 


44.73 


12.34 


66.19 


51.42 


RCC5 


34.14 


1.94 


53.95 


13.95 


98.10 


55.28 


RCU3 


34.33 


4.41 


76.02 


15.28 


102.81 


62.38 


RCU5 


36.78 


4.36 


78.14 


15.90 


112.62 


65.00 
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The first crack loads were obtained by visual examination. At this stage, the strengthened beams exhibit a 
ease of 174% compared to the control beams. The yield loads were obtained corresponding to the stage of 
loading beyond which the load-deflection response was not linear. At the yield load level, the GFRP strengthened beams 
showed an increase upto 166% compared to the control beams. The ultimate loads were obtained corresponding to the stage 
of loading beyond which the beam would not sustain additional deformation at the same load intensity. At the ultimate load 
level, the strengthened beams showed a maximum increase of 170% when compared with the control beams. From the 
experimental results, it can be observed that, at all load levels, a significant increase in strength was achieved by externally 
bonded GFRP laminates. This increase may be attributed to the increase in tensile cracking strength of Concrete due to 
confinement by the laminates. For the A series beams of steel ratio 0.419%, the ultimate load increased by 23.51% and 
58.81% for 3mm and 5mm thick CSMGFRP laminated beams. For beams strengthened with 3mm and 5mm thick 
UDCGFRP laminates, the ultimate load increased by 70.56% and 88.22%. The CSMGFRP strengthened HSC beams 
exhibited an increase in deflection which varied from 58.95% to 122.38% at ultimate load level. The UDCGFRP 
strengthened HSC beams exhibited an increase in deflection u hich varied from 153% to 171% at ultimate load level. 
For the B series beams of steel ratio 0.628%, the ultimate load increased by 13.65% and 18.19% for 3mm and 5mm thick 
CSMGFRP laminated beams. For beams strengthened with 3mm and 5mm thick UDCGFRP laminates, the ultimate load 
increased by 63.64% and 86.37%. The CSMGFRP strengthened HSC beams exhibited an increase in deflection which varied 
from 15.82% to 170.36% at ultimate load level. The UDCGFRP strengthened HSC beams exhibited an increase in deflection 
which varied from 189.98% to 211.59% at ultimate load level. For the C series beams of steel ratio 0.905%, the ultimate 
load increased by 12.49% and 66.72% for 3mm and 5mm thick CSMGFRP laminated beams. For beams strengthened with 
3mm and 5mm thick UDCGFRP laminates, the ultimate load increased by 74.72% and 91.40%. The CSMGFRP 
strengthened HSC beams exhibited an increase in deflection which varied from 144.27% to 162.61% at ultimate load level. 
The UDCGFRP strengthened HSC beams exhibited an increase in deflection which varied from 196.34% to 208.78% at 
ultimate load level. 



4.2 Ductility of beams 



TABLE 5 Ductility Indices of Tested Beams 



RCU3 
RCU5 



3.72 
4.17 



■1.08 
4.09 



Energy 

4.16 



.32 



Ductility is considered as an important factor in designing of structures especiall} in the seismic prone areas. The 
ductility of a beam can be defined as its ability to sustain inelastic deformation without loss in load carrying capacity, prior 
to failure. The ductility values for the beams were calculated based on deflection and energy absorption. The deflection 
ductility values were calculated as the ratio between the deflection at ultimate point to the deflection at yield point. The 
energy ductility values were calculated as the ratio of the cumul.m i nerg; ibsorption it ultimati tagi t< the cumulative 
energy absorption at yield. The ductility indices for the tested beams are presented in Table 5. The deflection ductility for the 
strengthened beams showed a maximum increase of 94.36%. 
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4.3 Comparison of Experimental Results with FEM Results 

The load - deflection response for tested beams is presented in Figs. 1 8-20. The general behaviour of the finite 
element models show good agreement with observations and data from the experimental tests. The failure mechanism of a 
reinforced concrete beam is modeled quite well using FEA and the failure load predicted is very close to the failure load 
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Figure 19. Load - Deflection Response of 'B ' Series Beams 
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V. CONCLUSIONS 

n the experimental results the follow ing conclusions are drawn: 
Strengthening of HSC beams using GFRP laminates resulted in higher load carrying capacity. The percentage 
increase in ultimate load varied from 23.51% to 88.22% for GFRP strengthened HSC beams. 
The percentage increase in deflection at ultimate stage varied from 15.82 % to 211.59% for HSC beams 
strengthened with GFRP laminates. 

The GFRP strengthened HSC beams show enhanced ductility. The increase in deflection ductility varied from 
17.36% to 100.38%. 

GFRP strengthened beams failed in flcxural mode only. 

The general behaviour of the finite element models show good agreement with observations and data from the 
experimental tests. The failure mechanism of a reinforced concrete beam is modeled quite well using FEA and the 
failure load predicted is vcr_\ close to the failure load measured during experimental testing. 
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